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ABSTRACT 

Turbulence  data  obtained  over  natural  water  waves  were  analyzed  using 
joint  probability  distribution  and  conditional  means  methods.  These  data 
represented  conditions  when  the  waves  were  decaying  and  when  the  waves 
were  building.  In  both  cases,  significant  wave-induced  fluctuations  were 
identified  in  the  airflow.  All  features  of  the  velocity  fluctuations 
were  examined  for  two  levels  above  mean  water  level.  In  the  case  of  a 
decaying  wave  field,  decelerations  in  the  airflow  can  be  associated  with 
an  assumed  propagating  pressure  maximum  over  the  crest  of  the  wave.  Other 
than  this  deceleration,  the  airflow  appears  to  reflect  simple  streamline 
bending  over  the  mobile  irregular  wave  surface.  In  the  case  of  a  building 
wave  field,  velocity  fluctuations  appear  to  agree  with  those  predicted  by 
linear  wave  generating  theories. 
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I.   INTEODUCTION 

This  study  examines  turbulence  data  obtained  over  natural  water  waves 
in  order  to  determine  the  influence  of  the  waves  on  the  overlying  air. 
The  data  examined  are  simultaneous  measurements  of  velocity  and  temperature 
fluctuations  at  two  fixed  levels  above  mean  water  level. 

The  study  extends  a  previous  analysis  on  data  for  which  significant 
evidence  of  wave-induced  motion  in  the  airflow  was  obtained.  Joint  proba- 
bility density  and  conditional  mean  analysis  are  the  statistical  procedures 
used  to  extend  the  previous  interpretations. 

The  lack  of  good  descriptions  of  the  interaction  between  the  ocean 
surface  and  the  adjacent  air  layer  is  a  major  obstacle  preventing  solution 
of  a  wide  class  of  applied  geophysical  problems  such  as  weather  prediction 
and  general  circulation  studies.  A  suitable  description  requires  the 
application,  to  observational  data,  of  specific  statistical  methods  which 
could  reveal  properties  of  wave-induced  organized  motion. 

Current  problems  requiring  additional  knowledge  on  the  dynamics  of  the 
surface  layer  are  the  determination  of: 

1.  kinetic  energy  dissipated  by  boundary  layer  turbulence; 

2.  thermal  energy  gained  by  turbulent  transport  of  sensible  and 
latent  heat; 

3.  energy  transformed  by  vertical  velocities  due  to  the  horizontal 
distribution  (convergence)  of  surface  friction;  and 

4.  surface  winds  used  in  semi -empirical  wave  forecasting  models. 
Although  numerous  observational  investigations  have  been  conducted  to 

describe  the  structure  of  the  near  surface  layer,  the  majority  of  these 
measurements  have  been  over  land.  Therefore,  parameterization  of  the 
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boundary  layer  over  the  sea  has  been  based  on  overland  results.  It  is 
clear,  however,  that  a  mobile  irregular  wave  surface  receiving  kinetic 
energy  from  the  airflow  can  influence  the  turbulent  regime  within  the 
adjacent  airflow.  Therefore,  empirical  formulae,  valid  or  at  least 
tested  over  land,  may  not  be  valid  over  water. 

Examples  of  the  conflicting  results  obtained  from  strictly  empirical 
approaches  appear  in  Figure  1  in  which  various  proposed  drag  coefficient 
relations  are  presented  for  over  the  water.  Being  able  to  add  some 
dynamical  description  to  this  region  should  help  in  narrowing  ones1  choices 
for  this  empirical  coefficient. 

Data  used  in  this  study  were  previously  analyzed  by  Davidson  (1970) . 
However,  those  analyses  consisted  of  spectral  analyses  and  the  limitations 
imposed  on  the  interpretation  of  spectral  results  hindered  the  description 
of  several  important  properties  of  the  observed  organized  motion.  Proce- 
dures employed  in  the  present  study  are  well  suited  for  examining  organized 
motion  within  a  turbulent  regime.  The  statistical  procedures  (Holland, 
1968)  were  designed  to  reveal  systematic  structure  in  various  types  of 
meteorological  data  but  at  the  present  time  have  only  been  applied  to 
larger  scale  fluctuations.  This  study  is  the  first  attempt  to  apply  these 
procedures  to  small  scale  fluctuations. 


fi  I  x"0 


FIGURE  1.  Sunmaries  of  Drag  Coefficients  Versus  Wind  Speed 
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II.   THEORETICAL  CONSIDERATIONS 

A.  WIND  WAVE-COUPLING 

In  a  survey  of  the  understanding  of  the  airflow  over  water,  Stewart 
(1967)  forewarned  that  the  air-water  boundary  is  not  "just  another 
boundary  layer",  but  is  unique  in  that  the  interface  is  a  part  of  the 
solution  for  the  boundary  layer  airflow.  The  interface  is  free  to  move 
and,  as  observed,  it  takes  advantage  of  this  freedom.  Stewart,  evaluating 
existing  wave  generation  models  (Miles  and  Phillips) ,  indicated  that  in 
1967  there  were  still  insufficient  overwater  measurements  to  either  verify 
or  invalidate  existing  analytical  theories  of  the  interaction  between  a 
turbulent  shear  flow  and  the  underlying  water  surface.  He  reported,  how- 
ever, that  even  though  he  examined  a  very  large  number  of  overwater 
velocity  spectra,  in  no  case  were  there  sharply  defined  spectral  peaks 
corresponding  to  the  wave  spectra  peaks.  Other  studies,  however,  have 
revealed  the  existence  of  wave-induced  motion  in  the  air. 

Yefimov  and  Pososhkov  (1969)  reported: 

"Laboratory  experiments  indicate  that  surface  waves  induce  in 
initially  quiet  air  not  only  periodic  velocity  fluctuations  but  also 
stationary  flow  in  the  direction  of  wave  propagation.  The  velocity 
of  stationary  flow  is  considerably  higher  than  would  follow  if  we 
assumed  it  to  resemble  Stokes  wave  flow  in  the  water. " 

Volkov  (1969)  reported,  based  on  interpretations  of  spectral  results: 

"...  the  wave  surface  has  a  considerable  influence  on  all  of  the 
statistical  characteristics  of  turbulence  in  the  atmospheric  boundary 
layer  -  the  fluctuation  spectra  and  dispersion,  the  turbulent  fluxes, 
etc.  We  therefore  succeed  in  finding  some  typical  features  and 
regularities  of  this  influence  by  the  complex  method  of  investigating 
the  processes  of  interaction  between  the  atmospheric  boundary  layer 
and  the  agitated  sea  surface . " 
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Davidson  (1970),  also  interpreting  spectral  results,  concluded: 

"...  analyses  of  the  boundary  layer  over  a  majority  of  the  earth's 
surface  cannot  be  based  on  the  assumption  that  the  surface  is  a  rigid 
boundary.  The  influence  of  the  waves  on  the  adjacent  airflow  has  been 
found  to  be  too  significant  to  rely  on  over-land  data  to  describe  the 
boundary  layer  over  waves.  Reliable  empirical  relations  giving  the 
flux  of  momentum  as  a  function  of  the  wind  profile  and  density  gradi- 
ents will  have  to  be  formulated  from  data  obtained  over  the  ocean 
surface  itself." 

An  analytical  model  describing  coupling  between  a  wave  surface  and 
the  overlying  air  flow  was  first  proposed  by  Miles  (1957) .  In  this  model, 
the  wave-induced  perturbation  is  considered  in  a  coordinate  system  which 
is  stationary  with  respect  to  the  wave  (i.e.,  surf -rider's  coordinate 
system  moving  at  the  phase  speed  c) .  Streamlines  in  such  a  coordinate 
system  would  appear  as  in  Figure  2.  The  Miles'  critical  height  (dashed 
line)  occurs  at  the  level  where  the  wind  speed  is  equal  to  the  wave  speed. 
The  wave,  Figure  2,  is  stationary  in  this  coordinate  system;  therefore, 
below  the  critical  level  the  air  appears  to  move  to  the  left  since  it 
moves  slower  than  the  waves.  Above  the  critical  level,  where  the  wind 
speed  is  greater  than  the  phase  speed,  the  air  moves  to  the  right. 

The  Miles  mechanism  can  be  described  "as  a  positive  feedback  mechanism 
and,  as  such,  it  depends  on  the  pre-existence  of  the  wave.  Forces  within 
the  Miles  theory  leading  to  wave  growth  are  depicted  in  Figure  2. 
Region  A  is  a  region  of  negative  vortex  acceleration  and  Region  B  is  one 
of  positive  vortex  acceleration,  but  the  absolute  value  of  the  vortex 
acceleration  is  greater  at  point  A  than  at  B.  Therefore,  there  is  a 
negative  vortex  acceleration  at  that  level  hence  a  deceleration  of  the 
mean  wind  when  averaged  over  a  wave  length.  This  deceleration  results  in 
a  loss  of  momentum  from  the  airflow  which  contributes  to  the  momentum 
increase  of  the  wave. 
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FIGURE  2.   Streamlines  of  sinusoidal  perturbation  in  a 
shear  flow  near  Z=  Z  ;  (a)  kinematic  conse- 
quence, (b)  consequence  of  dynamics  of  critical 
level. 
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It  has  been  shown  (Kinsman,  1965,  p.  571)  that  the  kinematics  at  the 
critical  level,  described  above,  would  result  in  a  phase  shift,  with 
height,  of  the  wave-induced  vertical  velocity  component. 

Therefore,  it  is  observed,  within  a  linear  model,  that  the  wave- 
induced  fluctuations  can  be  important  in  altering  momentum  transfer  in 
the  air  layer  adjacent  to  the  waves.  Also,  specific  predictions  of  this 
theory  can  be  examined,  such  as  a  phase  shift  of  the  vertical  velocity 
with  respect  to  height.  It  should  be  noted,  however,  that  within  this 
investigation  possible  non-linear  interactions  between  the  wave-induced 
motion  and  the  shear  flow  will  also  be  examined. 

B.   STATISTICAL  CONSIDERATIONS  AND  PROCEDURES 

Statistical  procedures  used  in  this  study  were  described  by  Holland 
(1368).  Applied  to  three  variables,  these  procedures  yield  the  joint 
probability  density  function  (JDF)  for  a  pair  of  the  variables  with  a 
conditional  mean  function  (CMF)  of  a  third  variable. 

Holland  argued  convincingly  on  the  suitability  of  this  type  of 

analysis  for  turbulent  regimes.  In  reference  to  the  application  of  these 

procedures  to  a  turbulent  regime,  Holland  cited  Batchelor  (1953) : 

"It  is  a  premise  of  probability  theory  that  a  random  function  f  (a) , 
say,  defined  for  all  values  of  a,  is  determined  statistically  by  the 
complete  system  of  joint-probability  distributions  of  values  of  the 
function  at  any  n  values  of  a,  where  n  may  take  any  integral  value. 
Likewise,  the  infinite  field  of  turbulent  motion  is  determined 
statistically  by  the  complete  system  of  joint-probability  distribu- 
tions of  the  values  of  the  vector  velocity  u(x,t)  at  any  n  points  of 
of  space-time." 

A  joint  probability  density  function  is  a  means  to  represent  simul- 
taneously observed  values  of  pairs  of  variables  by  contours  on  a  two- 
dimensional  array.  A  trivariate  statistical  relationship  is  obtained  by 
computing  a  conditional  mean  of  a  third  variable  as  a  function  of  the  two 
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independent  variables  contributing  to  joint  probability  density  function 
(JDF) .  This  trivariate  statistical  relationship,  hereafter  denoted 
JDF/CMF,  can  therefore,  be  represented  by  two  sets  of  contours  on  a 
two-dimensional  array. 

Joint  probability  density  function  and  conditional  mean  function 
analysis  results  are  illustrated  in  Figure  3.  The  array  is  centered  on 
the  mean  value,  which  is  zero,  of  the  two  variables.  Each  cell  in  the 
row  or  column  consists  of  two  numbers.  The  bottom  number  represents  the 
JDF,  the  probability  (xlO3)  of  the  joint  occurrence  of  the  representative 
deviation  for  the  pair  of  variables,  and  the  upper  number  represents  the 
CMF,  the  mean  of  the  third  variable  corresponding  to  the  given  joint 
occurrence  of  the  pair. 

These  results  have  been  smoothed  by  a  nine-point  smoother  to  remove 
the  perturbations  viewed  as  being  physically  and  statistically  insignificant. 
This  smoothing  reduced  an  18  by  18  array  to  the  observed  16  by  16  array. 

The  following  list  of  steps  within  the  JDF/CMF  statistical  procedures 
will  be  described  in  the  ensuing  paragraphs: 

1.  normalization  and  computation  of  preliminary  statistics  of  the 
variables ; 

2.  bandpass  filtering; 

3.  computation  of  JDF/CMF; 

4.  checks  for  filter  and  random  noise;  and 

5.  computation  of  phase-amplitude  information. 

Computation  and  plotting  of  all  statistics  were  done  in  the  W.  R.  Church 
Computer  Facility1  of  the  Naval  Postgraguate  School. 


1The  analysis  utilized  the  IBM  360  Computer  and  a  CALCOMP  plotting 
device  for  the  computations  and  contour  plotting. 
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FIGURE  3.  JDF/CMF  Smoothed  Output  Array 
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1.  Normalization  and  Canputation  of  Statistics 

JDF/CMF  statistics  were  ccrnputed  from  normalized  values  of  the 
original  variables.  For  example,  if  x  (i)  =  values  in  one  of  the  three 
sets,  x  =  mean  of  all  x(i) ,  and  <x  =  standard  deviation  of  all  x(i) ,  the 
JDF  and  CMF  were  computed  for  (x(i)-x)/cr.  Normalized  values  were  obtained 
for  each  of  the  three  variables  included  in  a  JDF/CMF  analysis. 

In  addition  to  the  mean  and  standard  deviation,  other  statistics 
were  computed  and  tabulated  for  each  variable  and  pairs  of  variables  used 
in  the  analysis.  These  statistics  were: 

1.  variance; 

2.  skewness  of  each  variable; 

3.  kurtosis  of  each  variable; 

4 .  covariances ; 

5.  correlation  coefficients;  and 

6.  number  of  data  points  used. 

2.  Bandpass  Filter 

A  bandpass  filter  was  applied  in  order  to  isolate  the  significant 
frequency  band  (0.1Hz  to  0.8Hz)  as  indicated  by  previous  spectral  analysis 
on  the  data  (Davidson,  1970) .  There  were  several  numerical  filters  that 
could  have  been  used.  Investigation  revealed  that  the  most  desirable 
filter,  with  respect  to  sharpness  of  cutoff  and  magnitude  of  oscillations 
beyond  the  terminal  frequency,  was  an  inverse  transform  filter. 

The  response  curve  for  the  inverse  transform  bandpass  filter  with 
59  weights  is  shown  in  Figure  4.  The  filter  is  superimposed  (Figure  14) 
on  the  wave  spectra,  also  shown  are  the  velocity  variance  and  covariance 
spectra.  The  filter  passes  a  range  of  frequencies  from  about  0.1  Hz  to 
0.8  Hz. 
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The  application  of  this  bandpass  filter  enabled  analyses  of  wave- 
induced  fluctuations  without  contamination  by  higher  and  lower  frequency 
fluctuations . 


1.0- 


n  (Hz)  1 


FIGURE  4.  Response  Curve  for  Inverse  Transform  Filter,  59  Pt. 


3.  Computation  of  JDF/CMF 

A  JDF  represents  the  probability  of  a  joint  occurrence  for  a  pair 

of  variables.  Therefore,  for  each  pair  of  variables  two  indices  were 

computed  representing  the  position  of  simultaneously  observed  values  of 

the  variables  in  an  18  by  18  array  of  a/2  by  a/2  joint  class  interval. 

For  example,  if  the  first  values  in  the  two  sets  were  1.25  and  0.9 

standard  deviations  from  their  respective  means,  their  joint  occurrence 

would  be  tabulated  in  the  a/2  by  a/2  cell  with  indices  (3,2) .  The  CMF 

was  computed  for  each  joint  class  interval  and  was  the  average  value  of 

the  third  variable  in  the  set  for  all  observation  times  when  the  two 

independent  variables  had  values  within  the  joint  interval. 
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An  additional  step  in  these  procedures  was  to  divide  the  number 
of  joint  occurrences  by  the  total  number  of  observations.  This  gave  the 
probability  per  cr2/4  within  each  cell  of  the  18  by  18  array  which  was 
then  smoothed  to  a  16  by  16  array.  Likewise,  the  sums  of  the  third 
variable  were  divided  by  the  number  of  joint  occurrences  in  each  cell. 
This  yielded  the  average  deviation  from  the  mean  of  the  third  variable 
corresponding  to  a  particular  probability  per  a2   for  the  other  two 
variables . 

The  output  was  printed  out  in  a  16  by  16  array,  Figure  3,  and  the 
array  was  analyzed  by  Contour,  a  computer  library  subroutine,  and  plotted 
on  the  CALCOMP  plotting  device. 

An  analyzed  JDF/CMF  (for  a  random  number  field)  is  illustrated  in 
Figure  5.  The  contours  connect  equal  values  of  occurrences  in  the  JDF. 
They  are  concentric  circles  because  they  were  generated  from  a  random 
number  field  and  circles  indicate  (as  expected)  no  relationship  between 
the  two  variables  making  up  the  JDF.   A  statistical  relationship  between 
two  variables,  other  than  between  orthogonal  relations  such  as  a  sine  and 
cosine,  would  cause  asymmetrical  or  elliptical  contours.  The  heavier 
lines  in  this  two-dimensional  picture  represent  isolines  of  the  CMF  with 
the  heaviest  (dashed)  being  the  zero  line.  An  analysis  with  this  degree 
of  complexity,  containing  essential  relations  which  exist  among  three 
variables,  is  not  difficult  to  read  with  comprehension. 
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FIGURE  5.  Analyzed  JDF/CMF  For  Random  Numbers 

Any  pair  of  variables  can  be  used  to  compute  the  JDF  and  any 
variable  can  be  used  to  compute  the  CMF.  In  Figure  25,  the  JDF  is  u  and 
w  denoted  by  (u,w)  and  the  CMF  is  r\   as  a  function  of  u  and  w.  The  JDF/ 
CMF  is  denoted  by  n  (u,w) . 

4.  Checks  for  Filter  and  Random  Noises 

Bandpass  filtered  data  were  analyzed  and  compared  to  the  unfiltered 
data  to  insure  that  the  data  were  not  altered  by  the  bandpass  filter.  The 
analysis  of  the  same  variables  for  the  bandpassed  (A)  and  the  non-bandpassed 
(B)  data  are  shown  in  Figure  6.  It  can  be  seen  that  the  JDF  (asymmetrical 
circles)  for  the  two  sets  of  data  are  essentially  the  same.  CMF  features 
in  the  bandpassed  analysis  are  smoothed  versions  of  those  in  the  non- 
bandpassed  analysis.  The  resultant  differences  between  A  and  B  are  differ- 
ences that  could  be  caused  by  the  addition  of  random  numbers  (which  is 
essentially  the  contribution  from  background  turbulence)  onto  the  data. 
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W1.5 


W1.5 


FIGURE  6.  JDF/CMF  of  Filtered  (A)  Data  and  Unfiltered  (B)  Data 

As  a  test  on  the  results,  random  numbers  were  superimposed  on  the 
data  used  in  each  analysis.  The  values  of  these  random  numbers  were  of 
the  same  order  of  magnitude  as  the  data  field  itself.  It  was  expected 
that  by  superimposing  random  numbers,  as  in  Figure  5,  onto  a  data  field, 
such  as  Figure  17A,  where  there  is  some  statistical  relationship,  there 
would  result  a  final  JDF  like  Figure  17B  where  the  asymmetry  would  be 
decreased.  This  was  done  for  each  variable  analyzed  and  such  a  result 
was  indeed  realized. 

The  analysis  of  each  set  of  variables  was  done  twice,  once  with 
the  bandpassed  data  and  once  with  a  random  number  field  superimposed  on 
the  bandpassed  data.  This  allowed  a  comparison  between  the  data  and  the 
data  plus  random  numbers  to  examine  if  the  significant  modes  within  the 
results  represented  something  other  than  chance  occurrences. 
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5.  Theory  and  Computation  of  Phase  Amplitude 

In  this  study,  information  was  also  desired  on  the  statistical 
dependence  of  the  various  turbulent  variables  on  the  phase  of  the  waves. 
Within  the  context  of  a  JDF/CMF  analysis,  this  would  be,  for  each  variable 
plotted  (Figure  18) ,  its  conditional  mean  function  with  respect  to  the 
phase  angle  of  the  wave.  Similar  to  the  JDF/CMF  results  a  method  of  dis- 
play is  required  to  emphasize  particular  variations  "and  eliminate  unwanted 
background  noise. 

An  approach  suggested  by  Holland  was  used  in  which  the  polar 
coordinates  representing  the  amplitudes  and  phase  angles  of  a  variable 
can  be  determined  by  computing  the  JDF  for  a  variable  and  its  time  deriva- 
tive. This  method  for  computing  phase  relations  should  reveal  (Holland, 
1968)  if  the  phase  relation  between  tvvo  variables  shifts  as  a  function  of 
the  amplitude  of  the  reference  variable.  In  particular,  it  will  reveal  if 
the  amplitude  of  the  dependent  variable  varies  non-linearly  with  the  ampli- 
tude of  the  independent  variable.  Such  a  result  is  highly  significant  in 
this  study. 

Holland  defined  for  any  fluctuating  meteorological  variable  y  a 
"statistical  eddy  period",  x  ,  defined  as 

a 

t  =  2it— 
o    a. 

If  the  period  of  oscillation  is  x  ,  the  ratio  of  the  standard  deviations 
can  be  found.  If  the  variable  is  written  as 

y  =  A  sin  2tt— 


then 


2irA     n  t 
y  =  cos  2tt— 
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The  means  are  both  zero  and  the  standard  deviations  are 


A 

a     =  y= 


and 


a. 


2ttA 

y  Vt  i 


therefore,   it  can  be  seen  that 


a.        2ir 

y 

The  statistical  eddy  period  for  each  variable,  the  standard  devi- 
ation of  each  variable  and  the  standard  deviation  of  the  tine  derivative 
of  each  variable  are  given  in  Table  1. 


VAR. 

19   AUGUST 

26  SEPTEMBER 

^y 

*y 

r. 

ff» 

y 

r0 

n 

25.03 

8.10 

3.8  8 

8.91 

4.11 

2.7  2 

uL 

26.97 

10.02 

3.38 

39.92 

14.65 

3.42 

wL 

2  2.13 

7.91 

3.52 

17.78 

7.72 

2.89 

Tt 

0.04 

0.01 

3.22 

0.16 

0.06 

3.42 

Uu 

27.86 

9.3  0 

3.76 

41.56 

14.48 

3.61 

wu 

23.26 

8.67 

3.37 

26.68 

10.54 

3.18 

Tu 

0.04 

0.01 

3.13 

0.10 

0.04 

3.36 

-UWL 

636.25 

304.66 

2.62 

1720.55 

906.79 

2.38 

-uwu 

871.74 

363.76 

3.01 

1759.83 

820.39 

2.70 

WTL 

0.85 

0.42 

2.56 

3.23 

1.65 

2.4  7 

WTU 

0.87 

0.44 

2.50 

3.9  3 

1.83 

2.70 

TABLE  1.  Statistical  Eddy  Period  of  Each  Variable  and  Standard 
Deviations  of  Each  Variable  and  its  Time  Derivative 
for  19  August  and  26  September  1968 
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Phase  results  consisted  of  the  statistical  dependence  of  various 
turbulent  variables  on  the  phase  of  the  waves  (n) .  Therefore,  a  JDF 
coordinate  system  was  used  with  n  as  the  ordinate  and  its  derivative  (n) 
as  the  abscissa.  In  this  coordinate  system,  with  the  passage  of  a  wave, 
n  increases  upward  and  successive  values  in  time  are  higher  on  the  right 
and  lower  on  the  left  of  the  ri-axis. 

If  the  surface  wave  were  a  pure  sinusoid,  successive  values  of 
JDF  (n,n)  would  follow  a  counterclockwise  circular  path  around  the  JDF 
coordinate  system,  passing  through  the  extremes  of  r\  when  crossing  the 
n-axis  and  through  extremes  in  n  when  crossing  the  n-axis.  However,  water 
waves  are  not  pure  sinusoids  so  these  paths  would  not  necessarily  be 
circular.  They  would,  however,  be  less  complicated  than  passages  of 
eddies  in  a  random  turbulent  regime. 

The  phase  angle  was  defined  by  the  angle  measured  (by  octants) 
counter  clockwise,  Figure  7,  from  the  positive  f|-axis  in  the  (n,n) 
coordinate  system. 

Phase  relations  were  examined  after  separating  the  results  into 
wave  amplitude  classes.  The  classes  were  defined  by  a  regions  on  the  JDF 
array.  In  this  case  two  amplitudes  were  considered,  these  being  small 
amplitudes  (less  than  1.5a)  and  large  amplitudes  (greater  than  1.5a).  The 
numbered  octants  are  shown  in  Figure  7  in  which  the  box  around  the  center 
represents  1.5a.  Essentially  this  permitted  an  examination  of  phase 
relationship  with  respect  to  large  and  small  waves.  In  contrast,  spectral 
analysis  yields  phase  information  as  a  function  of  frequency  only,  the 
amplitude  being  the  mean  amplitude  at  that  frequency. 
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FIGURE  7.  JDF/CMF  Smoothed  Output  Array  Depicting  the  Octants  for 
Phase-amplitude  Computations  and  1.5a  Boundaries  (Box 
Around  Center)  for  Separating  Large  and  Small  Amplitude 
Waves . 
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Applying  concepts  of  the  above  discussion,  computations  of  the 
phase  amplitude  relations  were  made  using  the  dependent  variables  as  the 
conditional  mean  and  the  waves  and  time  derivative  of  the  waves  as  the 
independent  variables  for  the  JDF.  After  obtaining  the  JDF/CMF  smooth 
array  output,  for  example  Figure  3,  the  products  of  the  dependent  variable 
(the  upper  number  in  each  cell)  and  the  probability  function  (the  lower 
number  in  each  cell)  in  all  joint  intervals  of  each  octant  and  amplitude 
class  were  computed  to  obtain  the  amplitude  stratified  phase  conditional 
mean  function.  This  was  done,  Figure  18  bottom  to  top,  for  n#  u,  w,  T, 
-uw,  and  wT,  shown  as  horizontal  bars  for  each  octant. 

A  smooth  curve  was  drawn  through  each  set  of  horizontal  bars  so 
as  to  make  the  two  areas  between  each  bar  and  the  curve  approximately 
equal.  These  smooth  curves  were  then  analyzed  for  their  phase  relation- 
ship with  the  wave  curve  (bottom  curve) .  The  left  hand  set  of  curves  (A) 
represents  the  small  amplitudes.  The  middle  set  (B)  represents  the  large 
amplitude  and  the  right  hand  set  (C)  is  the  total  of  A  and  B. 
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III.   DESCRIPTION  OF  DATA 

The  data  examined  in  this  study  were  frcm  two  days  (19  August  and 
26  September  1968)  of  measurements  on  Lake  Michigan  by  Davidson  (1970) . 
These  were  simultaneous  measurements  from  a  fixed  tower,  of  wind  (u  and 
w)  and  temperature  fluctuations  and  waves.  The  velocity  and  temperature 
measurements  were  obtained  at  two  different  levels  above  the  mean  water 
level.  Waves  were  measured  at  a  point  directly  below  the  air  sensors. 

A.   LOCATION  AND  DESCRIPTION  OF  THE  TOWER2 

The  measurements  were  made  on  Lake  Michigan  from  a  U.S.  Lake  Survey 
Research  tower  located  in  15  meters  of  water  1.6  kilometers  from  shore 
near  Muskegon,  Michigan  (Figure  8) .  Open  water  fetches  at  the  tower 
location  vary  from  80  to  300  miles  from  the  south  clockwise  through  the 
northwest.  Because  the  lake  is  not  so  large  that  separate  weather 
systems  and  hence  swell  propagate  from  different  directions,  the  local 
wave  field  was  related  to  the  wind  field  over  the  entire  lake. 

The  general  construction  of  the  research  tower  and  its  instrumentation 
is  shown  in  Figure  9.  The  structure,  an  open  triangular  truss,  was 
designed  both  to  achieve  stability  and  to  minimize  disturbance  to  the  air- 
flow and  water  motions.  These  were  important  considerations  in  the  present 
analysis  and  interpretation.  Stability  was  achieved  by  supporting  the 
tower  on  the  lake  bottom  by  a  tripod  which  was  secured  by  a  concrete  anchor 
at  each  leg.  The  tower  was  also  guyed  in  order  to  withstand  extreme  wind 
and  wave  conditions. 


2An  abridged  version  of  the  description  by  Davidson  1970) 
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B.   INSTRUMENTATION  .AND  SENSOR  MOUNTING  ARRANGEMENT 

Figure  10A  shows  the  features  of  the  sensor  mounting  arrangement.  A 
plumb  line  arrangement  of  sensors  enabled  determination  of  phase  relations 
between  the  waves  and  turbulent  components  in  the  air. 

The  turbulent  velocity  measurements  were  made  with  a  constant- 
temperature  hot-wire  anemometer  using  two  wires  in  an  X-configuration. 
The  two  wires  arranged  in  X-configuration  allowed  the  system  to  simul- 
taneously measure  the  two  wind  components  of  interest  (u  and  w) .  The  probes 
were  calibrated  both  in  situ,  based  on  simultaneous  cup  and  propeller 
anemometer  measurements,  and  in  a  wind  tunnel. 

A  resistance  thermometer,  in  a  Wheatstone  bridge  circuit,  was  used  to 
measure  the  temperature  fluctuations.  The  output  signal  of  the  system  was 
125  millivolts  for  a  one  degree  (Celcius)  temperature  change,  Testing 
this  system  in  a  laboratory  indicated  it  had  the  resolution  to  measure 
temperature  to  an  accuracy  of  0.05  C. 

The  waves  were  measured  by  a  capacitance  gauge  and  bridge  system, 
Figure  11.  The  bridge  was  adjusted  for  full  deflection,  zero  to  one 
milliamp,  between  no  immersion  and  full  immersion. 

The  velocity  and  temperature  sensors  and  wave  gauge  extended 
westward  2.5  times  the  width  of  the  tower's  side.  Extension  arms  consis- 
ted of  five  centimeter  rectangular  tubing  and  were  sturdy  enough  to  avoid 
motion  due  to  the  airflow  or  waves.  Also,  a  vertical  pipe  extended  from 
below  the  water  to  two  meters  above  the  water  so  the  sensors  could  be 
positioned  at  any  height  from  the  wave  crest  to  two  meters.  Above  the 
two  meter  level,  arms  were  located  at  the  four,  eight,  and  fifteen  meter 
levels . 
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Figure  10B  is  a  photograph  of  the  sensors  mounted  on  the  vertical 
pipe  above  the  mean  water  level  (MWL) .  A  wind  vane,  such  as  that  shown, 
was  used  at  all  levels  to  monitor  the  direction  of  the  mean  wind  with 
respect  to  the  sensors.  The  wave  gauge  can  be  seen  directly  below  the 
hot  wires  in  Figure  10B. 
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FIGURE  11.   Picture  of  capacitance  wave  gauge  and  bridge' 
amplifier  system. 
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IV.   PRESENTATION  AND  INTERPRETATION  OF  RESULTS 

In  this  section,  JDF/CMF  and  phase  amplitude  results  are  presented 
and  interpreted  for  the  most  probable  relations  between  fluctuations  in 
the  airflow  and  the  waves.  Temperature  relationships  were  also  computed 
and  a  cursory  examination  has  revealed  some  interesting  fluctuations  in 
the  temperature  field,  but  a  detailed  analysis  of  the  temperature  relation- 
ships was  not  performed  in  this  study.  They  are  included  for  reference 
in  Appendix  B. 

Non-linear  properties  and  phase  relations  of  velocity  fluctuations 
are  emphasized  in  the  following  interpretation.  Spectral  results 
(Davidson,  1970)  for  these  same  data  will  be  shown  to  complement  the 
discussion.  Spectral  results  are  used  to  indicate  the  intensity  of  the 
observed  wave-induced  fluctuations  with  respect  to  the  background,  "normal" 
turbulent  regime. 

Interpretations  obtained  from  the  analyses  of  the  present  study  should 
be  viewed  as  significant  extensions  to  those  available  from  the  previous 
spectral  results.  Although  the  spectral  results  show,  according  to 
frequency,  the  existence  of  wave  induced  fluctuations,  interpretations 
are  limited  with  respect  to  non-linear  properties  of  the  fluctuations  and 
phase  relations  between  fluctuating  variables.  Non-linear  properties  are 
not  available  within  spectral  results  because  of  the  inherent  linear 
expansion,  i.e.,  Fourier  Series.  Also,  phase  information,  obtained  from 
spectral  analysis,  is  highly  uncertain  because  of  the  high  level  of 
background,  random  turbulence  at  all  the  frequencies. 
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For  each  period,  two  figures  are  used  to  show  JDF/CMF  results  and 
two  figures  are  used  to  show  phase  amplitude  results.  Each  of  the  two 
figures  representing  JDF/G4F  results  include  a  trivariate  probability 
distribution  from  the  following  considerations  of  three  variables, 
(n  =  waves,  u  =  upper  level,  L  =  lower  level) , 

1.  nfo^) 

2.  n(uu,wu) 

within  each  of  these  figures,  six  plates  (consisting  of  three  pairs) 
are  used  to  show  the  three  possible  combinations  of  the  three  variables. 
The  pairs  arise  because  all  results  are  computed  with  and  without  the 
inclusion  of  random  numbers.  For  example,  six  plates  corresponding  to 
the  trivariate  (1)  above  (where  n  (u,w)  means  JDF  of  u  and  w  with  n  as 
CMF)  are: 

Without  Random  i-Tumbers  With  Random  Numbers 

1.  n  (u,w)  4.  n  (u,w) 

2.  w(n,u)  .  5.  w(ri/U) 

3.  u(n,w)  6.  u(n/W) 
Phase-amplitude  results  are  presented  in  two  figures,  one  figure  for 

all  variables  at  each  level  and  the  wave.  Each  of  these  figures  include 
the  following  arrangement  of  the  variables: 

1.  wT     Top  Panel  in  Figure  18 

2.  -uw 

3.  w 

4.  u 

5.  n     Bottom  Panel  in  Figure  18 
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A.  RESULTS  FOR  19  AUGUST 

This  period  was  from  1226  to  1244  CST.  Measurements  in  the  air  were 

made  at  1.5  and  4.0  meters  above  mean  water  level.  Higher  waves  existed 

in  this  period  than  during  the  other  period  considered,  26  September. 

The  following  descriptions  for  general  conditions  were  given  by  Davidson 

(1970) : 

"The  weather  map  for  0700  and  the  hourly  wind  and  wave  histories  for 
the  day  appear  in  Figure  12  (of  the  present  study) .  A  warm  front  had 
previously  passed  over  southern  Lake  Michigan  earlier  in  the  morning. 
As  a  result  of  this  frontal  passage  the  winds  were  southerly  and  the 
atmospheric  stratification  was  stable.  The  temperature  difference 
between  the  4.0  meter  level  and  the  surface  was  1.2  C  to  1.5  C.  From 
mid-morning  to  late  afternoon,  there  were  scattered  clouds.  Wind  (at 
the  16  meter  level)  and  wave  histories  (H1/3  cm,  hourly  averages) 
indicate  that  the  waves  began  to  develop  under  the  influence  of  the 
steady  south  wind  following  frontal  passage.  Even  though  the  wind 
increased  during  the  day,  the  wave  field  appears  to  have  obtained 
equilibrium  by  mid-moming  (0900)  and  remained  so  throughout  the 
observation  period.  The  wind  direction  and  wave  propagation  direction 
were  observed  to  coincide. 

Figure  13  shows  wind  (2  and  8  meters) ,  temperature  (degrees  Kelvin) , 
and  wave  (H1/3)  conditions  for  the  observational  period.  During  the 
period  there  was  a  slight  increase  in  wind  speed  but  the  wave  field 
appears  to  have  decreased  during  the  last  half  of  the  period.  The 
dynamics  of  the  matched  layer  should  have  had  no  influence  on  the 
results  because  the  critical  level  was  effectively  at  infinity." 

Variance  and  covariance  spectra  for  this  period,  obtained  by  Davidson 
(1970),  appear  in  Figure  14.  Velocity  variance  and  covariance  spectra 
for  the  1.5  meter  and  4.0  meter  levels  appear  on  sides  A  and  B  respectively, 
The  wave  spectrum,  along  with  the  bandpass  filter  (dashed  line) ,  appears 
on  both  sides  in  order  to  relate  the  wave  spectrum  peak  to  extrema  in  the 
velocity  spectra. 

Primary  features  in  these  spectra  are  extrema  in  both  variance  and 
covariance  spectra  near  0.2  Hz,  the  wave  spectrum  peak.  These  occur  for 
both  levels.  Energy  maxima  appear  near  this  frequency  in  the  variance 
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FIGURE  12.  General  Conditions  for  19  August  1968 
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WIND  AND  WAVE  HISTORIES 
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FIGURE  13.  General  Conditions  for  1226  to  1244  CST,  19  August  1968 
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FIGURE  14.  Velocity  Variance  and  Covariance  Results  and  Wave  Spectra 
(With  Bandpass  Filter  Superirtposed)  for  19  August  1968, 
With  1.5  Meter  Level  on  the  Left  and  4.0  Meter  Level  on 
the  Right 
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spectra  of  both  components  at  both  levels  near  or  centered  on  this 
frequency.  Covariance  spectrum  extrema  for  1.5  meters  indicate  an  upward 

transfer  of  momentum,  <j>   (n)>0,  and  that  for  4.0  meters  indicate  relatively 

uw  J 

larger  negative  values  for  <J>   (n)  f  resulting  in,  perhaps,  enhanced 
downward  transfer  of  momentum. 

Because  these  results  were  obtained  at  a  fixed  level  above  mean  water 
level,  energy  concentrations  at  the  frequency  of  the  wave  spectrum  peak 
could  reflect  simple  bending  of  streamlines,  as  depicted  in  Figure  ISA, 
in  a  coordinate  system  moving  at  the  speed  of  the  surface  wave.  In  this 
realization,  observed  fluctuations  in  u  and  w  would  occur  as  shown  in 
Figure  15B  and  correspond  to  those  observed  in  potential  flow.  Referring 
to  Figure  15B,  u  and  w  are  in  quadrature,  90°  out  of  phase,  and  would  not 
contribute  to  the  shear-induced  momentum  flux. 

in  uie  roxxowing  paragrapns ,  propercxes  or  uig  wave— induced  KtO"cicn  m 
the  airflow  will  be  examined  in  order  to  determine  why  contributions  to 
the  momentum  transfer  were  observed  and  why  opposite  contributions  were 
observed  at  two  levels  along  a  plumb  line.  Comparable  relations  for  a 
linear  process,  shown  in  Figure  15,  suggest  that  the  observed  spectral 
results  are  due  to  non-linear  properties  of  the  wave-induced  fluctuation 
within  the  shear  flow.  Dynamics  of  the  critical  level  were  not  a  factor 
in  the  results  for  this  period  because  the  critical  level  was  well  above 
the  level  of  measurements,  i.e.,  the  base  speed  of  0.2Hz  wave  is  equal  to 
7  m  sec-1  and  the  windspeed  at  4.0  meters  was  5  m  sec-1. 
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FIGURE  15.  Streamlines  for  Potential  Flow;  (A)  Simple  Bending  of 

Streamlines,  (B)  Relations  Between  u  and  w  in  Potential 
Flow 
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JDF/CMF  results  for  u,  w,  and  the  waves  (n)  appear  in  Figures  16  and 
17  for  the  1.5  and  4.0  meter  levels  respectively.  A  general  comparison 
between  plates  A  in  these  two  figures  reveals  a  significant  difference 
between  the  u-w  JDF  contours  at  the  two  levels.  U-w  JDF  contours  at  1.5 
meters,  Figure  16,  are  close  to  being  concentric  circles  while  those  for 
4.0  meters,  Figure  17,  are  asymmetrical  with  a  major  axis  135° 
counterclockwise  from  the  positive  u-axis. 

Although  results  appearing  in  plate  A  of  these  figures  represent, 
statistically,  the  probability  of  joint  occurrences  of  fluctuations  of 
various  magnitudes  of  u  and  w,  further  discussion  of  the  results  will  be 
within  the  context  of  physical  processes.  Therefore,  these  results  will 
be  interpreted  from  the  point  of  view  that  the  respective  quadrants, 
counterclockwise  from  the  positive  u-axis,  represent  the  following 
processes i 

1.  Quadrant  I  -  upward  transfer  of  positive  momentum; 

2.  Quadrant  II  -  upward  transfer  of  deficit  momentum; 

3.  Quadrant  III  -  downward  transfer  of  deficit  momentum;  and 

4.  Quadrant  IV  -  downward  transfer  of  positive  momentum. 
Herein,  negative  fluctuations  in  u  represent  deficit  momencum  and 

positive  fluctuations  represent  positive  momentum  relative  to  the  mean 
momentum  at  the  observation  level. 

General  patterns  occurring  in  plate  A  agree  with  those  physical  pro- 
cesses interpreted  from  the  cospectral  results.  The  1.5  meter  cospectrum, 
with  the  band  defined  by  the  bandpass  filter,  indicated  that,  in  the  mean, 
there  was  a  reduction  in  the  net  downward  transfer  of  momentum.  The  reason 
for  this  (Figure  16,  plate  A)  is  probably  because  there  is  a  decrease  in 
downward  transfer  of  positive  momentum  (Quadrant  IV)  ana  an  increase  in  the 


42 


a > 


1 — > 


FIGURE  16.  U-W  JDF/CMF  Results  for  1.5  Inters  for  19  August  1968 
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FIGURE  17.  U-W  JDF/CMF  Results  for  4.0  Meters  for  19  August  1968 
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downward  transfer  of  deficit  momentun  (Quadrant  III) .  The  cospectrum 
(Figure  14)  for  the  4.0  meter  level  represented  enhanced  net  downward 
transfer  of  momentum.  This  is  complemented  by  the  JDF  results  in  which 
the  probability  contours  are  asymmetrical,  elliptical,  and  represent 
significantly  greater  probability  for  upward  transfer  of  deficit  momentum 
(Quadrant  II)  and  downward  transfer  of  positive  momentum  (Quadrant  IV) . 

For  the  1.5  meter  level  slight  asymmetry  of  low  probability  contours 
are  indicative  of  greater  than  normal  occurrences  of  larger  negative 
fluctuations  in  u  at  this  level.  These  large  negative  fluctuations  in 
u  do  not  appear  to  have  a  preferred  association  with  either  positive  or 
negative  fluctuations  in  w.  The  latter  interpretation  follows  from  the 
fact  that  the  probability  contours  have  a  flat  edge  appearance  at  negative 
u  values. 

The  role  of  the  waves  in  producing  the  above  general  patterns  and 
deviations  can  be  examined  from  CMF  contours  for  the  wave  heights  in 
plates  A  of  Figures  16  and  17.  The  general  CMF  patterns  agree  with  those 
expected  for  measurements  at  a  fixed  level  within  bending  streamlines, 
potential  flow.  Observed  fluctuations  in  u  and  w  with  respect  to  the 
waves  for  potential  flow  appear  in  Figure  15. 

For  potential  flow,  the  CMF  patterns  for  the  waves  would  be  maximum 
positive  wave  values  associated  with  negative  u  and  zero  w.  Maximum 
negative  wave  values  would  be  associated  with  positive  u  and  zero  w.  A 
sufficient  deviation  from  the  pattern  is  the  bimodal,  or  split,  maxima 
and,  to  a  lesser  extent,  minima  in  the  wave  CMF  distribution.  There  are 
sufficient  points  contributing  to  the  split  maxima  to  maintain  this 
pattern  even  with  the  addition  of  random  numbers  (Figure  16B) . 
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The  statistical  significance  of  the  bimodal  extrema  can  also  be 
evaluated  by  examining  the  JDF's  for  u  and  the  waves  and  w  and  the  waves 
which  appear  in  plates  C  and  E  of  Figure  16.  The  asymmetry  expected  in 
potential  flow,  due  to  negative  correlation  between  u  and  n t   does  not 
appear  in  plate  C  of  this  figure.  There  exists  a  very  slight  asymmetry 
along  the  135°  axis  but  it  is  not  as  large,  for  example,  as  that  between 
uw  appearing  in  plate  A  of  Figure  16.  In  fact,  the  JDF  for  u  and  n  are 
nearly  as  concentric  as  those  for  w  and  n  r   which  in  potential  flow  should 
be  concentric  because  the  two  variables  are  90°  out  of  phase. 

Probable  causes  for  the  deviations  in  both  JDF  and  CMF  results,  from 
that  predicted  by  potential  flow  theory,  will  be  examined  from  phase- 
amplitude  results.  Phase-amplitude  relations  for  the  variables  u,  w,  T, 
-uw,  and  wT  with  respect  to  the  wave  appear  in  Figures  18  and  19  for  the 
1.5  and  4.0  meter  levels  respectively-  As  previously  indicated  each 
figure  has  three  panels  which  represent: 

1 .  small  waves  (A)  -  waves  less  than  1 .  5a ; 

2.  large  waves  (B)  -  waves  greater  than  1.5a;  and 

3.  all  waves  (C) . 

In  Figure  18  the  wave  appears  to  be  a  near  perfect  sine  function  but 
the  traces  for  u  and  w  are  distorted.  Distortion  in  u  and  w  are  most 
evident  within  the  negative  fluctuation  of  each.  In  general,  at  both 
levels,  extrema  in  w  occur  over  the  node  being  positive  over  the  rising 
part  of  the  wave  and  negative  over  the  sinking  part.  In  general,  extrema 
in  u  are  180°  out  of  phase  with  the  wave  being  negative  over  the  crest 
and  positive  over  the  trough. 
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FIGURE  18.  Phase-Anplitude  Results  for  1.5  Meters  for  19  August  1968 
(A)  Small  Waves,  (B)  Large  Waves,  (C)  All  Waves 
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FIGURE  19.  Phase-Amplitude  Results  for  4.0  Meters  for  19  August  1968 
(A)  Small  Waves,  (B)  Large  Waves,  (C)  All  Waves 
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A  significant  deviation  from  pure  potential  flow  is  the  distortion 
in  the  u  phase-amplitude  for  the  1.5  meter  level.  Negative  fluctuations 
in  u,  at  1.5  meters,  extend  further  behind  the  wave  than  predicted  by 
potential  flow  theory.  In  fact,  the  maximum  negative  value  of  u  occurs 
behind  the  crest. 

Distortions  which  appear  insignificant  within  the  phase-amplitude 
depiction  for  u  and  w  are  emphasized  in  the  -uw  phase-amplitude  results 
which  will  now  be  considered.  Figure  20  is  included  in  order  to  compare 
the  observed  results  with  those  predicted  by  a  linear  expansion,  i.e., 
Fourier  series,  of  the  fluctuations.  Phase-amplitudes  represented  in 
Figure  20  were  obtained  with  sine  waves  using  phase  angles,  with  respect 
to  the  wave,  obtained  from  the  spectral  analysis.  The  resulting  symmetry 
along  the  wave  is  evident  in  this  figure.  Also,  there  is  a  mean  downward 
niomentum  transport  (-uvr-o)  when  the  phase  relation  between  u  and  w  is 
greater  than  90°. 

In  comparison  to  the  linear  predictions  shown  in  Figure  20,  -uw  phase- 
amplitudes  obtained  in  this  study  were  asymmetrical.  The  -uw  phase- 
amplitude  results,  for  large  waves,  panels  B  in  Figures  18  and  19,  reveal 
that  the  possible  processes  (upward  and  downward  transfer  of  positive  and 
negative  momentum)  are  not  as  predicted  by  results  from  a  linear  description, 
In  fact,  downward  transfer  of  positive  momentum,  denoted  as  point  A,  is 
much  smaller  than  any  of  the  other  processes. 

In  summary,  these  results,  which  correspond  to  a  case  where  the  waves 
were  in  a  state  of  equilibrium  or  slight  decay,  indicate  that,  as  a  first 
approximation,  the  airflow  adjacent  to  the  waves  conforms  to  the  irregular 
mobile  wave  surface.  This  was  particularly  true  for  the  1.5  meter  level 
where  JDF  contours  were  nearly  concentric.  However,  more  detailed  analyses 
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Analysis  Phase  Angles;  (A)  1.5  Meters,  (B)  4.0  Meters  for 
19  August  1968 
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of  JDF/CMF  results  and  also  analyses  of  phase -amplitude  results  revealed 
significant  distortions/  from  the  potential  flow  prediction,  in  the  u 
component.  This  was  not  the  case  for  the  w  component  which  exhibited, 
at  both  levels,  phase  relations  predicted  by  potential  flow. 

Phase-amplitude  results  for  -uw  appeared  to  show  most  clearly  the 
processes  associated  with  a  wave-induced  fluctuation  in  the  airflow.  At 
both  levels  upward  and  downward  transfer  of  deficit  momentum  appeared 
to  be  the  dominant  process.  Therefore,  assuming  that  random  or  shear- 
induced  properties  are  averaged  out  during  phase-amplitude  computations, 
these  results  indicate  that  the  waves  influence  did  not  enhance  the 
downward  transfer  of  positive  momentum. 

A  probable  cause  for  the  distortion  in  the  u  component,  at  both  levels, 
appears  in  Figure  21.  It  is  suggested  that  a  high  pressure  region, 
centered  slightly  behind  the  crest  and  propagating  at  the  speed  of  the 
wave,  produced  the  observed  decelerations,  and  hence  distortions,  from 
the  potential  flow  predictions  in  the  u  component.  Such  a  pressure 
pattern  coincides  with  the  general  condition  of  a  decaying  wave  field. 
Because  the  pressure  extrema  are  moving,  it  is  apparent  that  this  influence 
on  the  wind  field  is  not  the  same  as  that  by  a  stationary  object.  These 
interpretations  suggest  that,  for  this  case,  the  drag  coefficient  should 
approximate  the  smooth  flow  case  with  some  adjustment  for  the  observed 
pressure  influence. 

B.  RESULTS  FOR  26  SEPTEMBER 

This  period  was  from  1355  to  1413  CST  and  the  measurements  were  at 
1.5  and  4.0  meters.  Although  the  waves  were  not  as  high  during  this 
observation  period  as  they  were  on  19  August,  they  increased  slightly 
during  the  period.  The  wave  spectrum  was  observed  to  have  changed  sha'pe 
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FIGURE  21.  Summary  Illustrating  Probable  Cause  for  the  Distortion 
in  the  u  Component  for  19  August  1968 
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during  the  day.  This  change  was  a  shift  of  the  spectral  peak  to  a 
lower  frequency.  The  weather  for  the  day  (abridged  from  Davidson,  1970) 
is  shown  in  Figure  22.  The  wind  during  the  observation  period  was 
relatively  constant  and  the  hydrostatic  conditions  were  unstable  with 
the  water  temperature  approximately  5  C  higher  than  the  air  temperature 
at  four  meters.  Wind  and  wave  histories  are  shown  in  Figure  23. 

Variance  and  covariance  spectra  for  this  period  obtained  by  Davidson, 
appear  in  Figure  24  in  which  the  velocity  spectra  for  the  1.5  meter  and 
4.0  meter  levels  appear  on  sides  A  and  B  respectively.  Again  the  wave 
spectrum  peak,  with  the  bandpass  filter  (dashed  line) ,  appears  on  both 
sides . 

Similar  to  19  August,  the  primary  feature  in  these  spectra  are 
extrema  in  both  variance  and  covariance  spectra,  at  both  levels,  near 
the  wave  spectrum  peak  approximately  0.26  Hz.  Covariance  spectra  indi- 
cate enhanced  downward  transfer,  4   (n)<0,  of  momentum  for  both  levels 

uw 

contrary  to  the  19  August  spectra  which  showed  an  upward  transfer  of 
momentum  at  the  1.5  meter  level  at  0.2  Hz. 

JDF/CMF  results  for  u,  w,  and  n  appear  in  Figures  25  and  26  for  the 
1.5  meter  and  4.0  meter  levels  respectively.  A  general  comparison  of 
the  two  levels  (plate  A)  reveals  that  they  are  similar.  That  is,  the 
u-w  JDF  in  both  have  a  major  axis  of  asymmetry  approximately  135°  counter- 
clockwise from  the  positive  u-axis.  This  general  pattern  agrees  with  the 
physical  processes  interpreted  from  the  cospectral  results  which  suggest 
a  net  downward  transfer  of  momentum. 
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FIGURE  25.  U-W  JDF/CMF  Results  for  1.5  Meters  for  26  September  1968 
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FIGURE  26.  U-W  JDF/CMF  Results  for  4.0  Meters  for  26  September  1968 
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It  would  be  difficult  to  differentiate  the  JDF  results  from  those 
for  a  "normal"  shear  flow.  In  a  shear  flow,  u-w  JDF  probability  contours 
would  also  be  asymmetrical,  elliptical,  and  would  represent  significantly 
greater  probability  for  upward  transfer  of  deficit  momentum  (Quadrant  II) 
and  downward  transfer  of  positive  momentum  (Quadrant  IV)  .  An  apparent 
deviation  is  that  Quadrant  II  and  IV  are  not  symmetric.  Upward  transfer 
of  deficit  momentum  (Quadrant  II)  is  the  dominant  of  the  two.  Wave  CMF 
features  and  phase-amplitude  relations  will  be  the  primary  acts  in 
identifying  the  influence  of  the  waves.  At  the  1.5  meter  level,  joint 
occurrences  between  large  negative  fluctuations  in  u  and  positive  fluctu- 
ations in  w  are  strongly  associated  with  the  wave  crest.  This  relation 
is  not  as  evident  at  the  4.0  meter  level.  The  CMF  patterns  at  both 
levels  are,  of.  course,  not  the  same  as  those  predicted  by  potential  flow 
theory  (described  in  the  19  August  discussion) .  Similar  to  19  August, 
panels  C  and  E  do  not  show  significant  asymmetries  between  either  u  or  w 
and  the  waves.  However,  the  lower  level,  1.5  meters,  does  have  the  higher 
correlations . 

The  manner  of  the  deviations  from  potential  flow  is  quite  apparent  in 
the  phase -amplitude  results,  Figures  27  and  28.  Similar  to  19  August, 
the  wave  appears  to  be  a  nearly  perfect  sine  function.  At  both  levels,  u 
appears  to  have  the  same  general  relation  as  in  potential  flow,  negative 
over  the  crest  and  positive  over  the  trough.  This  relationship  is  almost 
perfect  at  1.5  meters.  The  4.0  meter  level  is  distorted  during  the 
positive  cycle  with  the  maximum  u  shifted  forward  of  the  trough  and  has 
a  flat  appearance.  Although,  the  phase  of  extrema  in  u,  for  4.0  meters, 
is  discernible  in  these  results,  the  amplitude  is  very  rmall  and,  in  fact, 
barely  above  a  level  which  would  be  viewed  as  significant. 
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FIGURE  27.  Phase-Amplitude  Results  for  1.5  Meters  for  26  September  1968 
(A)  Small  Waves,  (B)  Large  Waves,  (C)  All  Waves 
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FIGURE  28.  Phase-Amplitude  Results  for  4.0  Meters  for  26  September  1968 
(A)  Small  Waves,  (B)  Large  Waves,  (C)  All  Waves 
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The  primary  phase-amplitude  result  is  the  shift  in  the  w  maximum. 
Whereas  in  potential  flow,  and  on  19  August,  the  w  maximum  is  ahead  of 
the  wave  crest  (zero  over  the  crest) ,  in  this  observation  period  w 
maximum  is  directly  over  the  wave  crest. 

The  deviations  are  again  emphasized  in  the  -uw  trace  and  Figure  29 
is  included  for  comparative  purposes.  Figure  29  depicts  the  product  of 
two  sine  waves  whose  phase  angle,  with  respect  to  the  wave,  are  the 
same  as  u  and  w  and  which  were  obtained  from  spectral  analysis.  It  can 
be  seen  by  comparing  Figure  28  and  29  that  the  momentum  transfer  processes, 
as  depicted  by  JDF/CMF  results,  differ  significantly  from  the  results, 
depicted  by  a  linear  description.  The  dominate  process  (large  scale 
waves  (B)  in  Figures  27  and  28)  appears  to  be  upward  transfer  of  deficit 
momentum  and  the  only  other  process  of  any  significance  appears  to  be 
downward  transfer  of  positive  momentum.  Downward  transfer  of  positive 
momentum  is,  however,  much  less  than  the  upward  transfer  of  deficit 
momentum. 

In  summary,  these  results,  representing  a  case  in  which  the  dynamics 
of  the  critical  level  and  wave  generation  processes  could  be  a  factor, 
reveal  significant  differences  from  those  observed  in  the  other  period 
examined.  U-W  JDF  results  for  both  levels  indicated  that  upward  transfer 
of  deficit  momentum  and  downward  transfer  of  positive  momentum  were  the 
most  significant  processes.  These  same  results  indicated  that  the  former 
was  the  most  dominant  of  these  two  processes. 

Phase-amplitude  results  indicated  that  the  asymmetrical,  shear  flow 
type,  pattern  of  the  u-w  JDF  were  due  to  a  near  180°  phase  difference 
between  u  and  w.  This  phase  difference  existed  because  u  maintained  out- 
of  phase  relations  with  the  wave  extrema  but  the  maxima,  a  in  w  shifted, 
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FIGURE  29.  Phase-Amplitude  Results  for  Sine  Waves  Utilizing  Spectral 
Analysis  Phase  Angles;  (A)  1.5  Meters,  (B)  4.0  Meters  for 
26  September  1968 
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relative  to  the  19  August  description,  to  a  position  over  the  crest. 
The  shift  in  phase  of  the  w  component  is  predicted  in  the  linear  theory 
and  is  due  to  the  dynamics  at  the  critical  level. 

The  prominence,  again,  of  deficit  momentum  transfer  was  evident  in 
the  -uw  phase  relation.  A  physical  interpretation  for  this,  as  was 
possible  for  19  August,  has  not  been  achieved.  However,  it  could  be 
related  to  the  slight  distortions  in  both  the  u  and  w  phase-amplitude 
traces.  It  is  known  that  distortions  of  this  type  are  associated,  in  a 
linear  expansion,  with  super-position  of  different  harmonics.  It  is, 
therefore,  suggested  that  in  this  case  the  distortions  were  due  to  a 
coupling  between  the  wave-induced  fluctuations  and  a  fluctuation  of  a 
higher  frequency.  The  latter  could  be  due  to  buoyant  processes  which 
existed  over  this  relatively  warm  undulating  wave  surface. 
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V.  CONCLUSIONS 

In  this  study,  turbulence  data  obtained  over  waves  were  examined  to 
obtain  some  understanding  of  properties  and  influence  of  wave-induced 
fluctuations  on  the  turbulent  regime.  Statistical  procedures  were  used 
which  enabled  interpretations  to  be  made  on  non-linear  properties  of 
the  wave-induced  fluctuations  and  on  separate  modes  of  momentum  transfer. 
Two  data  periods,  representing  decaying  and  generating  wave  conditions, 
were  examined  and  the  results  revealed  significant  difference  in  these  two 
periods.  The  usefulness  of  the  selected  statistical  procedures  is 
supported  by  the  following  conclusions  from  the  results  of  each  of  these 
two  periods. 

For  a  case  when  the  near  surface  wind  speed  is  much  less  than  the 
speed  of  the  predominant  surface  wave,  the  critical  level  nearly  at 
infinity,  these  results  indicate: 

1.  The  primary  feature  of  the  flow  is  one  of  simple  streamline 
bending  in  which  the  airflow  conforms  to  the  wave  surface  and  this  can 
be  observed  as  high  as  4.0  meters  above  mean  water  level. 

2.  There  is  a  slight  deceleration  in  the  airflow  due  to  a  pressure 
maximum  located  near  the  wave  crest.  Because  this  pressure  maximum  is 
traveling  at  the  speed  of  the  wave  it  does  not  represent  an  obstacle  to 
airflow  as  a  stationary  object  would. 

3.  An  observed  increase  in  the  negative  correlation  between  hori- 
zontal and  vertical  components,  which  is  normally  considered  to  represent 
enhanced  downward  transport  of  positive  momentum,  was  due  to  enhanced 
upward  transfer  of  deficit  momentum. 
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4.  If  these  results  are  representative  of  conditions  when  the  waves 
are  in  equilibrium  with  the  wind  field,  an  appropriate  drag  coefficient 
would  be  one  reflecting  smooth  flow  conditions.  This  is  based  on  the 
fact  that  the  airflow  simply  follows  the  wave. 

For  a  case  where  the  near  surface  wind  speed  and  the  phase  speed  of 
the  predominant  surface  waves  are  about  the  same,  the  critical  level  near 
the  level  of  the  observations,  the  results  indicate  that: 

1.  There  is  an  increased  negative  correlation  between  u  and  w  that 
is  due  to  a  90°  phase  shift  of  the  wave-induced  w  component.  This  shift 
is  predicted  by  a  linear  theory  which  is  based  on  the  dynamics  at  the 
critical  level. 

2.  There  appears  to  be  a  distortion,  from  potential  flow  predictions, 
in  the  u  component  and  this  could  be  due  to  wave-induced  buoyant  motion. 

Future  research  on  these  data  and  other  data  obtained  over  natural 
water  waves  should  include  provisions  for  the  statistical  procedure 
employed  in  this  study.  Wave  following  sensors  (for  example  Boston, 
et.al.,  1969,  Davis,  1969)  should  be  used  in  the  measurements.  This 
would  remove  the  necessary  comparisons  with  potential  flow  predictions. 
Relations  between  temperature  fluctuations  and  the  waves  should  be 
interpreted. 
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APPENDIX  A:      JDF/CMF  COMPUTER  PROGRAM 
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/urcuNUiA  ±5:      TEMPERATURE  JDF/CMF  ANALYSES 
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FIGURE  30.  JDF/CMF  Results;  W1  5  and  T±       with  Waves  f  r  19  August  1968. 
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FIGURE  31.  JDF/CMF  Results;  W    and  T    with  Waves  .or  19  August  1968 
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FIGURE  32.  JDF/CMF  Results;  U,  c  and  U,  n  with  Waves  for  19  August  1968. 
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FIGURE  33.  JDF/CMF  Results;  W,  c  and  W,  „  with  Waves  or  19  August 
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FIGURE  34.     JDF/CMF  Results;   T,    _  and  T,   A  with  Waves  for  19  Aucrust  1968. 
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FIGURE  35.  JDF/CMF  Results;  W,   and  T,  5  with  Waves  for  26  Septerrter  1968, 
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FIGURE  36.  JDF/CMF  Results;  W.  Q  and  T.  Q  with  V?aves  for  26  September  1968. 
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FIGURE  37.  JDF/CMF  Results;  U,  5  and  U.  Q  with  Waves  for  26  September  1968, 
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FIGURE  38.  JDF/CMF  Results;  W,  _  and  W.  Q  with  Waves  for  26  September  1968. 
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FIGURE  39.  JDF/CMF  Results;  T.  _  and  T.  Q  with  Waves  f<  :  26  September  19C3. 
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